joined = geopandas.sjoin(stations, districts[['district_name', 'geometry']], op='within')